arXiv:astro-ph/0112042vl 3 Dec 2001 


Spatially resolved imaging spectroscopy of T Tauri 

M. E. Kasper 1 , M. Feldt, T. M. Herbst, S. Hippier 
Max-Planck-Institut fur Astronomie, Konigstuhl 17, 69117 Heidelberg, Germany 

mkasperOeso.org 
and 

T. Ott, L. E. Tacconi-Garman 

Max-Planck-Institut fur extraterrestrische Physik, Garching, Germany 

ABSTRACT 

The small separation of the individual stellar components in the TTau system has 
prevented spectroscopy of the individual components in the near-infrared (NIR) in the 
past. We present NIR data (H- and K-band) on this object taken using the ALFA 
adaptive optics (AO) system in combination with the 3D integral held spectrograph. 
Except for the Brackett series seen in emission, the NIR spectra of TTauS appear 
featureless, suggesting that warm dust dominates the radiation. No ro-vibrational lines 
of molecular hydrogen were found in close vicinity to the stars. The Brackett line 
emission from TTauN and S could not be resolved with the 3.5-m-telescope, implying 
that it arises within 6 AU from the stars. The line ratios between Bry (2.166 fim) and 
BrlO (1.736^m) toward TTauN and TTauS are similar, suggesting similar selective 
extinctions towards the respective emitting regions. Our results are consistent with a 
model that describes T Tau S as a pre-main sequence star surrounded by a small edge-on 
disk, leaving the polar regions relatively unobscured. We present numerical simulations 
which support this model. 

Subject headings: instrumentation: adaptive optics, integral held spectroscopy — bina¬ 
ries: pre-main sequence — stars: individual (T Tau) — infrared: stars 


1. Introduction 

T Tauri is a complex system of multiple (at least three) stellar components, jets, and outflows 
unlike that seen in any other T Tauri star. 
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The optically visible primary, TTauN, is an early K star (Cohen & Kuhi 1979), which has 
a companion at 0.7" separation to the south found by speckle interferometry in the NIR (Dyck 
et al. 1982). This infrared companion, TTauS, has not yet been detected at visible wavelengths 
(Stapelfeldt et al. 1998). TTauS exhibits a spectral energy distribution (SED) suggestive of a very 
young embedded source with a luminosity of about 11 Lq, making it somewhat more luminous than 
TTauN (Koresko et al. 1997). It shows an apparent lack of cold dust emission, implying a rather 
small mass for a circumstellar disk (Hogerheijde et al. 1997). TTauS underwent a brightness flare 
of two magnitudes at various wavelengths between 1989 and 1991 (Ghez et al. 1991; Kobayashi et 
al. 1994). After a quiescent low brightness phase in the mid 90s (Herbst et al. 1997), another flare 
started in 1998 (Roddier et al. 2000) and lasts until today. 

By speckle holography observations from 1998, Koresko (2000) found a companion to T Tau S 
at a separation of only 0.05" and assigned the name TTauSb to this component. Only two years 
later, this companion was found at a distance of 0.08" with its position angle changed by 28° with 
respect to TTauS (Kohler et al. 2000). Orbital motion studies of TTauN and S using speckle 
interferometry (Ghez et al. 1995) and AO (Roddier et al. 2000) imply a minimum mass of 3.3 M 0 
for the whole T Tauri system. 

The near environment of T Tauri is a source of surprisingly strong extended ro-vibrational 
emission of molecular hydrogen (Beckwith et al. 1978). At higher resolution, the spatial distribution 
of this excited H 2 shows two separate outflows, perpendicular to each other in projection, and rather 
strong emission in close vicinity to both stars (Herbst et al. 1996). The excitation appears to be 
produced by shock heating rather than UV fluorescence. It could either be caused by infall of 
material onto the star, or by outflows intercepting ambient matter. The Bry emission appeared to 
be confined to TTauN. In 1996/1997, one year after Herbst’s observations, TTauS showed weak 
Bry emission, while H 2 was barely detectable (Beck et al. 2001). 

Recent optical spectroscopy (Solf & Bohm 1999) indicated that the outflows seen in H 2 are in 
fact almost perpendicular with the north-south component seemingly originating from T Tau S at 
79° inclination from the line of sight, and the east-west component originating from TTauN. The 
latter coincides with the HH 155 jet, which is seen at an inclination angle of 23° (Eisloffel & Mundt 
1998). 


2. Observations 

The observations described here were taken with the ALFA adaptive optics (AO) system 
(Kasper et al. 2000) on the 3.5-m telescope at the German-Spanish Astronomical Center on Calar 
Alto. Using TTauN (my = 9.6) as a natural guide star, ALFA delivered angular resolutions 
(FWHM) of 0.14" in H- and 0.16" in K-band throughout the observations. The integral field 
spectrometer 3D yields spectra of an entire 16 x 16 pixel field at resolutions of R ~ 1000 to 
R « 2100 (A/AA) (Weitzel et al. 1996; Anders et al. 1998). The image scale of 0.07" per pixel 



- 3 - 


results in a field of view (FOV) of one arcsecond. Due to the small FOV, sky frames have to be 
observed in addition, effectively doubling the observing time. 

T Tauri was observed with ALFA and 3D on September 28, 1999 in K-band and on October 
1, 1999 in H-band with spectral resolutions of R ~ 1100 and R ~ 2100 respectively. The total inte¬ 
gration time on source was 18 and 14 minutes, respectively. Standard techniques of sky subtraction 
and flat fielding removed background emission, dark current, and pixel-to-pixel gain variations. 
The G2V star HD 26710 was observed at a similar airmass to T Tauri and served as a photometric, 
spectroscopic, and PSF reference. Division by the spectroscopic standard and subsequent multipli¬ 
cation by the sun’s spectrum at the same spectral resolution eliminated telluric and photospheric 
features of the G2V star and produced the correct spectral slope across the band. The final 3D data 
structure consists of an image cube with 600 separate wavelength slices, each containing spatial 
information over the FOV. This data cube allows one either to create line maps by subtracting the 
continuum next to a line from the line flux itself, or to extract a spectrum anywhere in the field. 


3. Results 
3.1. Spectra 

Viewing the 3D data cube signal in a direction perpendicular to the image planes produces 
conventional spectra. The traces in Figure 1 represent the total flux in a Q.25" square box centered 
on T Tau N and T Tau S, respectively. Results for T Tau S include also its companion because it was 
not resolved by our observations. The vertical bars and labels denote diagnostic lines in the actual 
photometric band. The detection limit for the lines can be estimated by measuring equivalent 
widths (EWs) on the noisy continuum. In the case of our spectra, such EWs are well below 0.15 A, 
suggesting that line features with EWs greater than about 0.3 A are detectable at a 2cr level. 

Table 1 lists the line fluxes and uncertainties based on least-squares fits of Gaussian profiles to 
the lines drawn from the spectra. The model includes five parameters: the central wavelength Ao, 
the amplitude A , the standard deviation of the Gaussian <r, and two parameters m and b giving 
the slope and intercept of the linear baseline, respectively. With the model 

F\ = b + m\ + A exp - ^ ^ , (1) 

2cr- 

the flux is given by \p2/nA ct, and the formal uncertainty is A F\ = [27t(<t 2 AA 2 + A 2 Acr 2 )] 1//2 , where 
A A and Act are those changes in the A and a parameters which increase the y 2 per degree of 
freedom by one (Bevington & Robinson 1992). It should be noted that a number of additional 
effects which are not included in the uncertainties have an impact on the line measurements. These 
effects include, for example, imperfect telluric calibration or uncertainty in the brightness of the 
photometric standard. 

No significant H 2 emission could be found within one arcsecond of the stars. Brackett recorn- 
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bination lines of atomic hydrogen dominate the H- and K-band spectra of T Tau N and T Tau S 
with comparable EWs for both stars. Flux ratios between the Bry and BrlO line are about 3.25 
for both TTauN and TTauS, consistent with the value expected for optically thin emission from a 
gas in local thermal equilibrium and obeying Menzel Case B recombination (see Osterbrock 1989, 
p. 84). 

Photospheric absorption features as given in Table 2 are only found for T Tau N. None are 
detected in the spectrum of TTauS. The relative strengths of the EWs are consistent with those 
seen in early K spectral standard star spectra (e.g. Meyer et al. 1998; Wallace & Hinkle 1997), 
but the absolute values are much more shallow. This requires a large amount of continuum veiling 
such as the contribution of the hot inner part of a circumstellar disk. In the case of T Tau N, the 
K-band excess must be of the order 80 % of the total flux. There is some hint in both stars of 
the presence of the v = 2 — 0 vibrational bandhead of CO near 2.3 /ini. This wavelength region 
suffers from telluric absorption features, however, and is difficult to correct completely. Also, there 
is no evidence of additional Av = 2 bandheads (3 — 1, 4 — 2, etc.) at longer wavelengths, which 
accompany the 2 — 0 bandhead in spectra of late type stars (Kleinmann & Hall 1986). 


3.2. Linemaps 

The continuum subtracted line maps displayed in Figure 2 allow us to search for extended 
Brackett emission when compared to images in the adjacent continuum. The spatial resolution of 
the H-band line maps from October 1999 is about 0.14" FWHM. When compared to the presumably 
point-like continuum PSF, a potential source extension in a Brackett line image would add to this 
value. We performed empirical tests using Gaussian fits to the PSF cores in the continuum linemaps 
next to the BrlO line. We found that, at the distance of T Tau (~ 140 pc), a source extension of 
more than 6 AU convolved with the mean PSF is required to increase the FWHM significantly (3(7) 
compared to the scattering in the continuum PSF FWHMs. The continuum subtracted BrlO PSF 
does not show a significantly larger FWHM than the continuum PSFs, meaning that the radiation 
arises within 6AU from the stars. Using a similar test on the position showed that the Brackett 
emitting region is coincident with the star at a level of 0.8 AU. 


4. Discussion 

Apart from the Brackett emission lines, the H- and K-band spectra of TTauS are featureless. 
In comparison to our detection limit of about 0.3 A, typical EWs of atomic features for late type 
standard stars are of the order of a few Angstroms (Meyer et al. 1998; Wallace & Hinkle 1997). 
This implies that the stellar photosphere of T Tau S is heavily veiled, if visible at all, meaning that 
we basically observe hot dust. A scenario, where TTauS is a star reddened by large amounts of 
extinction from cold dust, can thus be ruled out. 
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Our data do not show any H 2 emission within one arcsecond of TTauN and S. This is a 
somewhat surprising result, because Herbst et al. (1996) found strong H 2 emission centered at the 
positions of TTauN and TTauS, and Beck et al. (2001) detected weak H 2 emission from TTauS 
only. Brackett recombination lines instead are prominent in the spectra of both stars, with the 
contribution of TTauS being stronger than observed by Herbst et al. (1996) who did not detect 
Bry emission and Beck et al. (2001) who found a lineflux of 1.6 ± 0.5 x 10^ 13 erg cm -2 s -1 . Our 
observations date from a flare period, while the two others were carried out during a low brightness 
phase of TTauS. Hence, the variable emission line luminosity of TTauS suggests a common cause 
with its variable brightness. 

In an attempt to estimate the excess luminosity produced by an accretion event from the Bry 
luminosity, we used the following empirical formula given by Muzerolle et al. (1998) 

log(L acc /L 0 ) = (1.26 ± 0.19) log(L Br7 /T 0 ) + (4.43 ± 0.79), (2) 

to derive a value of L acc = 0.44 ± O.O9L 0 for TTauS. The uncertainty includes our measurement 
error as well as those related to Equation 2. The derived value is certainly not sufficient to explain 
brightness flares of T Tau S of more than one magnitude when taking into account its total lumi¬ 
nosity of more than 10 L 0 (Koresko et al. 1997). The larger Bry flux during the brightness flare 
is however consistent with a foreground extinction that has decreased from Ay ~ 9 to Ay ~ 6 in 
the transition between quiescent and elevated flux state (Beck et al. 2001). These extinction val¬ 
ues were derived from observed water ice absorption and therefore represent foreground extinction 
produced by relatively cold material. They do not represent the effects of warm material close to 
the star. 

The extinction towards a Brackett emitting region can be determined by the observed Bry to 
BrlO flux ratio x, if the intrinsic recombination line flux ratio xo is known. Although our observed 
flux ratio of x = 3.23 is consistent with Menzel Case B, this could be a coincidence since Evans et 
al. (1987) did not find the corresponding hydrogen recombination line ratios in his sample of pre- 
main-sequence objects. Moreover, Natta et al. (1988) showed that relative intensities of hydrogen 
recombination lines created in a mostly neutral stellar wind, instead of a fully ionized region, are 
a complex function of wind parameters and stellar radiation field. Let us, however, suppose that 
the intrinsic recombination line ratios of TTauN and TTauS are similar, neglecting the actual 
generation process. Since the extinction towards T Tau N is negligible in the NIR (Cohen & Kuhi 
(1979), Ay = 1.39), an intrinsic xq = 3.23 is adopted. With the NIR extinction law from Mathis 
(1990), A(A) oc A~ L7 , it follows that A Br io/A Br7 = 1-46 or A Br io — A Br7 = 0.46A Br7 . Hence, the 
observed line ratio depends on extinction as 

x = x 0 ■ io-°- 4 (°- 46Ab ^). ( 3 ) 

The maximum flux ratio for T Tau S consistent with the accuracy of our observations (see Table 1) is 
x = 4.1, yielding an upper limit on the extinction towards the Brackett emitting region of Ay ~ 6. 
The star instead must be must be obscured by at least Ay = 7 (Stapelfeldt et al. 1998). The 
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conclusion is that in contrast to the star, the Bry emitting region around T Tau S does not seem to 
be heavily obscured. 

Although, we observed at the diffraction limit of the 3.5-rn telescope, the Brackett emitting 
region around TTauN and S could not be resolved. This is consistent with the neutral wind 
model of Natta et al. (1988) that predicts most of the hydrogen recombination line flux to be 
produced within a few tenths of an AU around the central star, a scale not resolved by our data. 
This result does not contradict radio observations showing TTauS as an extended source of two 
circularly polarized lobes approximately 0.15 arcseconds apart (Ray et al. 1997), because it is 
unlikely that this radio emission is produced by free-free transitions. Almost certainly, its origin is 
gyrosynchrotron radiation which would not be detectable in our Brackett linemaps. 


5. Model calculations 

A proto-stellar or Class I model with an infalling envelope seems unlikely for TTauS, mostly 
because of a lack of cold dust in its immediate environment and the presence of its close companion. 
Based on the orientation of T Tau S’s jet almost perpendicular to the line of sight and the proximity 
to TTauSb, Koresko (2000) proposed a scenario in which both stars are surrounded by small disks 
seen nearly edge-on. This could easily provide the required local extinctions without large amounts 
of dust and could further explain the brightness variations by changes in a disk which processes 
stellar photons via a combination of extinction and scattering. 

What can our data add to this picture? The lack of photospheric absorption features deduced 
from the H- and K-band spectra rules out a pure extinction scenario. The relatively low extinction 
towards the Brackett-series emitting region indicated by our data further suggests a geometry, 
where the disk is observed at an inclination such that it obscures the star while leaving the line of 
sight towards the polar regions comparatively unobscured. NIR flux could be contributed by the 
hot inner parts of the disk itself, if the optical depths in the disk are not too large. 


5.1. The model 

To test the validity of these assumptions, we performed model calculations using the radiative 
transfer (RT) code of Manske et al. (1998). This code is based on a method described in detail by 
Men’shchikov & Henning (1997). The model consists of a small disk with an opening angle of 30° 
inclined by 14°, i.e. almost edge-on. Hence, the line of sight towards the central star intersects the 
disk just beneath its surface. The model parameters are given in table 3. The code was used in its 
2D mode, calculating at first the temperatures and mean intensities in the disk’s midplane and at 
its conical surface. Intermediate temperatures and mean intensities are computed by interpolation. 
Finally, specific intensities are computed using a set of rays intersecting the disk. Based on these 
specific intensities, the SED is computed. 
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5.2. Result 

The SED of T Tau S as obtained from the model calculations is given in Figure 3. To account 
for the observed foreground extinction of about Ay = 6 mag during an elevated flux state (see Sec. 
4), the result of the RT computation is shown reddened by this amount using the reddening law 
from Mathis (1990). 

The modeled SED well reproduces the flux densities measured by Ghez et al. (1995) during 
the 1990 flare. Only the data point at 20 //,m could not be fitted by the small disk, because of a 
lack of dust at larger distances to the star and therefore lower temperatures. However, the missing 
20 /un-flux could be produced by the same dust that provides the foreground extinction. While 
little is known about geometry and temperature distribution of a circumbinary dust torus around 
TTauS, a rough estimate shows that such a torus located between 25 and 50 AU from the stars 
with a dust density of about 10 -14 Kgm -3 and a temperature of 75 K can produce the required 
foreground extinction and excess flux density of about 10 Jy at 20 /xm. 

The very small disk radius of 3 AU is on the one hand an input requirement due to the close 
companion TTauSb at a projected distance of about 11 AU. According to Artymowicz & Lubow 
(1994), a circumstellar disk in a binary system has a radius smaller than about 0.45 times the 
binary separation for a circular binary orbit. Larger eccentricities can reduce this value to 0.2 
times the semimajor axis of the binary. On the other hand, the small disk radius is required by the 
RT calculations, because larger disks tend to produce more MIR and less NIR flux when observed 
edge-on and therefore could not fit the data. 

In conclusion, the concept of a small truncated disk surrounding T Tau S with dimensions that 
fit the close binary system as a circumprimary disk and produce a high extinction along the line of 
sight, while leaving the Bry emitting region (at least the northern half of it) relatively unobscured, 
are not contradicted by the RT computations. Whether the details of the model are significant, 
and the above implications regarding the circunrsystem dust torus are correct, could be subject to 
further research. 


6. Summary 

Using AO compensated imaging spectroscopy, we were able to spatially resolve T Tau N and 
T Tau S and to obtain individual H- and K-band spectra of the two components for the first time. 
The observations were carried out during an elevated flux state of the variable TTauS. It was not 
possible to disentangle the contributions of T Tau S and T Tau Sb. 

H 2 emission was not detected in the immediate vicinity of TTauN and TTauS, whereas both 
exhibit the Brackett series in emission with comparable equivalent widths. Compared to previous 
observations of T Tau S during a quiescent phase, the Bry lineflux has increased significantly while 
H 2 emission could no longer be found. This is evidence for a common cause of variable emission 



line spectrum and luminosity. The accretion luminosities deduced from the Bry fluxes are typical 
for pre-main sequence stars and do not suggest exceptionally high accretion rates. 

The ratio between Bry and BrlO luminosity is similar for TTauN and TTauS. Considering 
the low extinction for TTauN and assuming a similar intrinsic ratio for both stars, the Brackett 
series emitting region around T Tau S does not appear as heavily obscured as the star. 

The Brackett emission lines of TTauN and TTauS arise within a few AU of the stars. They 
were not resolved by our observations. The apparent mismatch between Brackett- and radio 
linemaps are further evidence against a fully ionized regions as a source of the radio emission. 

Photospheric absorption features were only detected for T Tau N which appears heavily veiled. 
The absence of any photospheric features in the case of T Tau S rules out a pure extinction scenario 
as an explanation for its unusual SED. 

Given that we do not see the stellar photosphere of TTauS, but do see the Brackett series 
emitting region seemingly unobscured, a small edge-on disk may be consistent with the observed 
SED when viewed at an angle such that the star is just obscured. We carried out radiative transfer 
calculations using a very simple model and found that a small edge-on disk can indeed reproduce 
the SED of T Tau S. 


We are grateful to both the ALFA and 3D teams for their invaluable help, and the Calar Alto 
Observatory for their hospitality. We extend our thanks to Tracy Beck for providing a preprint 
of her paper, and to Volker Manske and Olaf Kessel for their support with the radiative transfer 
calculations. We also thank Andy Nelson for discussions on the mechanics of circumstellar disks in 
binary systems as well as editor and referee for their valuable comments to the manuscript. 
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Fig. 2.— Continuum subtracted line maps in BrlO and Bi’ 7 . The images display the whole FOV 
of 3D. A squareroot scaling was chosen to enhance the fainter TTauS. 
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Table 1: Emission line fluxes (in 10 13 erg cm 2 s 3 ) and EWs (in A) of TTauN and TTauS. 


Line 

Wavelength 

TTauN 

TTauS 


(m m ) 

Flux 

EW 

Flux 

EW 

H I 






Brl3 

1.611 

2.2 ± 1.0 

0.54 ± 0.25 

0.2 ±0.05 

0.45 ±0.1 

Brl2 

1.641 

0.6 ± 0.2 

0.16 ±0.03 

0.2 ±0.05 

0.45 ±0.1 

Bril 

1.681 

0.9 ± 0.6 

0.22 ±0.16 

0.4 ±0.08 

0.79 ±0.16 

BrlO 

1.736 

3.0 ± 0.6 

0.79 ±0.15 

0.8 ± 0.09 

1.42 ±0.16 

Br7 

2.166 

9.7 ± 1.1 

3.2 ±0.36 

2.6 ±0.29 

2.16 ±0.25 

h 2 






1-0 S(l) 

2.121 

a 





not detected 
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Table 2: Equivalent widths of atmospheric features of TTauN. 

Species Mg A1 OH Ca 

A 0 (pm) 1.575 1.711 2.282 1.674 1.688 2.263 

EW(A) 0.84 ±0.45 0.55 ± 0.21 0.33 ± 0.12 0.75 ± 0.44 0.51 ± 0.33 0.27 ±0.18 
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Table 3: Parameters used in the radiative transfer calculation to model the SED of TTauS. 


Central object 

Luminosity: 

16 L 0 


Temperature: 

5000 Kelvin 


Distance: 

140 pc 


Dust Distributions 

Density distr.: 

rsj y ^ 


Outer radius: 

3 AU 


Inner radius: 

0.9 AU 


Dust Properties 

Size distr.: 

a = 0.005 - 0.25 

/irn 


N(a) ~ a” 3 ' 5 


Silicate/Carbon ratio: 

Si:C 3:2 


Optical data : 

Preibisch et al. 

(1993) 


Draine & Lee (1984) 

Optical depth of disk: 

~700 at 0.5 n m 


Total dust mass: 

2.6 x 10“ 5 M o 


Foreground extinction: 

Ay = 6 mag 


Disk geometry 

Opening angle: 

30 deg 


Inclination: 

14 deg 













